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ABSTRACT: We investigate the two dimensional free 
convection and mass transfer flow of an incompressible, 
viscous and electrically conducting fluid past a continuously 
moving vertical flat plate in the presence of heat source, 
thermal diffusion, large suction and the influence of uniform 
magnetic field applied normal to the flow. Usual similarity 
transformations are introduced to solve the momentum, 
energy and concentration equations. To obtain the solutions 
of the problem, the ordinary differential equations are solved 
by using perturbation technique. The expressions for velocity 
field, temperature field, concentration field, skin friction, 
rate of heat and mass transfer have been obtained. The 
results are discussed in detailed with the help of graphs and 
tables to observe the effect of different parameters. 

Keywords'. Free convection, MHD, Mass transfer, Thermal 
diffusion, Heat source parameter. 

I. INTRODUCTION 

Magneto-hydrodynamic (MHD) is the branch of 
continuum mechanics which deals with the flow of 
electrically conducting fluids in electric and magnetic fields. 
Many natural phenomena and engineering problems are 
worth being subjected to an MHD analysis. Furthermore, 
Magneto-hydrodynamic (MHD) has attracted the attention 
of a large number of scholars due to its diverse applications. 

In engineering it finds its application in MHD 
pumps, MHD bearings etc. Workers likes Hossain and 
Mandal [1] have investigated the effects of magnetic field 
on natural convection flow past a vertical surface. Free 
convection flows are of great interest in a number of 
industrial applications such as fiber and granular insulation, 
geothermal systems etc. Buoyancy is also of importance in 
an environment where differences between land and air 
temperatures can give rise to complicated flow patterns. 

Convection in porous media has applications in 
geothermal energy recovery, oil extraction, thermal energy 
storage and flow through filtering devices. The phenomena 
of mass transfer are also very common in theory of stellar 
structure and observable effects are detectable, at least on 
the solar surface. The study of effects of magnetic field on 
free convection flow is important in liquid-metal, 
electrolytes and ionized gases. The thermal physics of 
hydro-magnetic problems with mass transfer is of interest in 
power engineering and metallurgy. The study of flows 
through porous media became of great interest due to its 
wide application in many scientific and engineering 
problems. Such type of flows can be observed in the 
movement of underground water resources, for filtration and 



water purification processes, the motion of natural gases and 
oil through oil reservoirs in petroleum engineering and so 
on. The porous medium is in fact a non -homogeneous 
medium. The velocity is usually so small and the flow 
passages are so narrow that laminar flow may be assumed 
without hesitation. Rigorous analysis of the flow is not 
possible because the shape of the individual flow passages is 
so varied and so complex. Poonia and Chaudhary [2] studied 
about the flows through porous media. Recently researchers 
like Alam and Rahman [3], Sharma and Singh [4] 
Chaudhary and Arpita [5] studied about MHD free 
convection heat and mass transfer in a vertical plate or 
sometimes oscillating plate. An analysis is performed to 
study the effect of thermal diffusing fluid past an infinite 
vertical porous plate with Ohmic dissipation by Reddy and 
Rao [6]. The combined effect of viscous dissipation, Joule 
heating, transpiration, heat source, thermal diffusion and 
Hall current on the hydro-Magnetic free convection and 
mass transfer flow of an electrically conducting, viscous, 
homogeneous, incompressible fluid past an infinite vertical 
porous plate are discussed by Singh et. al. [7]. Singh [8] has 
also studied the effects of mass transfer on MHD free 
convection flow of a viscous fluid through a vertical channel 
walls. An extensive contribution on heat and mass transfer 
flow has been made by Gebhart [9] to highlight the insight 
on the phenomena. Gebhart and Pera [10] studied heat and 
mass transfer flow under various flow situations. Therefore 
several authors, viz. Raptis and Soundalgekar [11], Agrawal 
et. al. [12], Jha and Singh [13], Jha and Prasad [14] have 
paid attention to the study of MHD free convection and 
mass transfer flows. Abdusattar [15] and Soundalgekar ET. 

al. [16] also analyzed about MHD free convection 
through an infinite vertical plate. Acharya et. al. [17] have 
presented an analysis to study MHD effects on free 
convection and mass transfer flow through a porous medium 
with constant auction and constant heat flux considering 
Eckert number as a small perturbation parameter. This is the 
extension of the work of Bejan and Khair [18] under the 
influence of magnetic field. Singh [19] has also studied 
effects of mass transfer on MHD free convection flow of a 
viscous fluid through a vertical channel using Laplace 
transform technique considering symmetrical heating and 
cooling of channel walls. A numerical solution of unsteady 
free convection and mass transfer flow is presented by Alam 
and Rahman [20] when a viscous, incompressible fluid 
flows along an infinite vertical porous plate embedded in a 
porous medium is considered. 
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In view of the application of heat source and 
thermal diffusion effect, the study of two dimensional MHD 
free convection and mass transfer flow past an infinite 
vertical porous plate taking into account the combined effect 
of heat source and thermal diffusion with a great agreement 
was conducted by Singh ET. al [21]. 

In our present work, the effect of large suction on 
MHD free convection heat and mass transfer flow past a 
vertical flat plate has been investigated by similarity 
solutions which are obtained by employing the perturbation 
technique. We have used MATHMATICA to draw graph 
and to find the numerical results of the equation. 

II. THE GOVERNING EQUATIONS 

Consider a two dimensional steady free convection 
heat and mass transfer flow of an incompressible, 
electrically conducting and viscous fluid past an electrically 
non-conducting continuously moving vertical flat plate. 

Introducing a Cartesian co-ordinate system, x-axis 
is chosen along the plate in the direction of flow and y-axis 
normal to it. A uniform magnetic field is applied normally to 
the flow region. The plate is maintained at a constant 
temperature T w and the concentration is maintained at a 
constant value C w . The temperature of ambient flow is 7^ 
and the concentration of uniform flow is Co. 

Considering the magnetic Reynold's number to be 
very small, the induced magnetic field is assumed to be 

negligible, so that B = (0, B 0 (x), 0).The equation of 

conservation of electric charge isV. / Where / = (Jx, Jy, 
Jz) and the direction of propagation is assumed along y-axis 
so that / does not have any variation along y-axis so that the 

dJ 

y derivative of / namely = 0 resulting in /^constant. 

dy 

Also the plate is electrically non-conducting J y =0 
everywhere in the flow. Considering the Joule heating and 
viscous dissipation terms to negligible and that the magnetic 
field is not enough to cause Joule heating, the term due to 
electrical dissipation is neglected in the energy equation. 

The density is considered a linear function of 
temperature and species concentration so that the usual 
Boussinesq's approximation is taken as p=p 0 (T-T^) 
+/?*(C-Co)}]. Within the frame work of delete such 
assumptions the equations of continuity, momentum, energy 
and concentration are following [21] 
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The boundary conditions relevant to the problem are 



--u n 



0 



(4) 



(5) 



= 0, 



= v 0 (*), T = T w C = C w at y- 
= 0, T = T a0t C = C ao at y^^ 

Where u and v are velocity components along x- 
axis and y-axis respectively, g is acceleration due to gravity, 
T is the temperature. K is thermal conductivity, o is the 
electrical conductivity, D M is the molecular diffusivity, U 0 is 
the uniform velocity, C is the concentration of species, B 0 (x) 
is the uniform magnetic field, C P is the specific heat at 
constant pressure, Q is the constant heat source(absorption 
type), D T is the thermal diffusivity, C(x) is variable 
concentration at the plate, v 0 (x) is the suction velocity, p is 
the density, 3 is the kinematic viscosity, ft is the volumetric 
coefficient of thermal expansion and ft* is the volumetric 
coefficient of thermal expansion with concentration and the 
other symbols have their usual meaning. For similarity 
solution, the plate concentration C(x) is considered to be 

C(jC) — Coo ( Cyv~ Cao)X. 
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III. MATHEMATICAL ANALYSIS suction 

We introduce the following local similarity X(£) = l + sX l (%) + s 2 X 2 (%) + s 3 X 3 (%) + ( 16 ) 

variables of equation (6) into equation (2) to (4) and 

boundary condition (5) and get the equation (7) to (9) and = sY i^) + £ Y 2 (£) + £ Y 3 (£) + 

boundary condition (10) z(£) = + s 2 Z 2 (%) + ^ 3 Z 3 (<J) + ( 18 ) 

Using equations (16)-(18) in equations (12)-(14) and 

y/ = ^23xU f(rj) rj = y I ^° ^) considering up to order 0(s 3 ), we get the following three 

V2i9x sets of ordinary differential equations and their 

Q( \ - ^ ~ ^ ' Ay( \ - ~ Q corresponding boundary conditions 

~ T w -T„' ~ C w -C„ First order 0(e): 

f'" + ff'-Mf f + Gr6> + Gm<p = 0 (7) Xf+Xf = 0 (19) 

0* + Pr/0'-5Pr0 = O (8) Y{+PiY{ = 0 (20) 

0" - 2Scf 0 + Scffl + S o Sc0" = 0 (9) Z{ + &Z1 + ScS 0 Y{= 0 (21) 

#/Q? The boundary conditions for 1 st order equations are 

where Pr = is the Prandtl number, ^ Tz 

ik ^=0,^ = 1, y i= l, Z 1= lat £ = 0 ( 22 ) 

Gr= 2xgfi(T„-TJ . s ^ Gmshofnumber9 ^ =1? y i=0? Zi=0 as ^ = 0 

^0 2 

Second order <9(e j: 

Gm= 2x8 ^ (C ^~ C " ) is Modified Grashof + Xj +A r 1 Z 1 " =MX{ -Gr^ -GmZ l (23) 

number, tf + ftXtf + MS = SPrY, (24) 

2xO + ScZ' 2 = 1ScX[Z x - ScK x Z{ - ScS 0 Y% (25) 

S = is the Heat source paramater, . 

U 0 The boundary conditions for 2 order equations are 

w 2xa% 2 (x) . „ . X 2 =0, X^=0, F 2 =0, Z 2 =0 at £ = 0 (26) 

M = is Magnetic parameter, 

U 0 p X 2 = 0, F 2 =0, Z 2 =0 as <J^oo 

e *9 . i n r • i, i , Third order Of/j: 

5c = is the Schmidt number and 

D M y X'"+ XI + X"X 2 X X X 2 = MX' 2 - GrY 2 - GmZ 2 (27) 

S Q =I^ is the Soret number *3 + PrX^ + PrX 2 Y{ + PrF 3 ' = SPrY 2 (28) 

c - c °° z; + scz 3 ' = 2Scx;z 2 + 2Scx 2 z 1 - scx 2 z[ (29) 

with the boundary condition, _ ^ _ y, t 

y* _ j f'=l 0 = 1 (j) = l at 77 = 0 (10) b° un( l ar y conditions for 3 rd order equations are 

T a ' , ' v X 3 =0, X 3 '=0, F 3 =0, Z 3 =0 at ^ = 0 (30) 

f =0, 0 = 0, 0 = 0 as 77^oc 3 3 3 

, X 3 '=0, F 3 =0, Z 3 =0 as ^ = oo 
| 2x 

where / w = v 0 (jc) and primes denotes the The solution of the above coupled equations (19) to (30) 

* 0 under the prescribed boundary conditions are as follows 

derivatives with respect to rj. First order 0(e): 

Z = r?f w , f(t]) = f w X(&, X x =\-eS ( 31 > 

0(r 1 ) = f w 1 Y^), 0(Tj) = fJZ(£) (ID Yy =e-^ ( 32 > 

By using the above equation (11) in the equations (7) - _ Sc ^ _ Pr ^ 

(9) with boundary condition (10) we get, Z 1 =(l-A 1 )e c + A x e 

X"' + XX" = s{MX' - GrY - GmZ) (12) Second order 0(s 2 ): 

1 (^4) 

Y" + PrXY' = sSPrY (13) X 2 = -e' 24 +(^ + A>£)*T # - \e** +V +Ao 

Z"-2ScZX' + ScXZ' + ScS 0 Y" = 0 (14) Y 2 =(A n - A^e'^ - A^ 1 ^ ^ 

The boundary conditions (10) reduce to „ , /0 ^. 

X-i, r-e, Y-s andZ.„ a, (-0 (15) 2. = (36) 

I'-o. r = o, z=o as ^ + 

Third order Ofe J: 

Now we can expand a, F and Z m powers or s as 

follows since g = is a very small quantity for large 

f 2 
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X 3 = B 31 + (5 24 + B^ + ^ <f )e~ f - - A, - *± 0 

e- 2 < + (B 25 + B 26 + B^)e- p " - [B 21 - B 28 - BJR 

1 



(37) 



-B 1Q e 



24 



B 



(38) 



-(1+Pr)£ 



+ B 37 e' 



-(2+Pr)| 



2 

Z 3 =B 122 e-^+( J B 116 +5 97 ^ 



+ B x e 



-(Pr+Sc)£ 



+ (^117 +^100^> 



-Q+Pr)£ 



-(2+Pr)£ 



(39) 



+ £ e~ {2+Sc)4 +B e~ (Fr+Sc) ^-B e~ 2V ^ 

-r ^ 103 e ^104^ ^105^ 



Using equations (16) to (18) in equation (11) with the 
help of equations (31) to (39) we have obtained the 
velocity, the temperature and concentration fields as 
follows 

Velocity distribution 

u = U 0 f'(? 7 ) = U 0 f 2 X'^) 

= u 0 [x;&+£X^)+s 2 x 3 ^)] 

Temperature distribution 

= Y l (Z) + £Y 2 ^) + s%(Z) 
Concentration distribution 
m = flz(g) = Z x {£) + sZ 2 (£) + ^ 2 Z 3 (£) 
The main quantities of physical interest are the local skin- 
friction, local Nusselt number and the local Sherwood 
number. The equation defining the wall skin-friction as 

(43) 



(40) 



(41) 



(42) 



y=0 



Thus from equation (40) the skin-friction may be written 

as 

* 00 [/"07)L> = [1 + £ V + A " 2 A + Apr 2 



+ A^Sc 2 } + 8 2 {B 9 - 2B 32 + B 24 
- B 7 + 2^ - 2£ 18 Pr+ Pr 2 (£ 25 
+ B 26 ) -2£ 21 5c - Sc 2 (B 21 - B 2% ) 

-£ 29 (1 + Pr) 2 -B 30 (l + Sc) 2 -l}] 



(44) 



The local Nusselt number denoted by Nu and defined as 

(45) 

Nu = - ~" 



Therefore using equation (41), we have 
M/ oc [0'(n)\^ = -Pr+ ^(A u - A 12 ) + ^ 2 [(Z? 41 
-PrZ? 40 -(l + Pr)(Z? 39 -Z? 35 ) 



(46) 



£ 36 -£ 37 (2 + Pr)-£ 38 (Pr 



+ 5c)-^Pr] 

The local Sherwood number denoted by Sh and defined 



as 



Sh- 



3C_ 



(47) 



(48) 



Hence from equation (42), we can write 
Sh oc [<f>'(7i)l =Q = - A Pr- Sc(l - A ) + s{Sc(B, + B 2 
+ B 3 -B 5 )-B 6 -Vr(B 3 -B 5 ) + B A 
-B^l + V^-B^l + Sc) 
+ £ 2 {-£ 122 Sc-(l + Sc)£ 116 + £ 97 
-(1 + Pr)5 117 +5 100 -5 102 (2 + Pr) 
-5 103 (2 + 5c)-Z? 104 (Pr+5c) 
+ 25 105 Pr+25 106 5c-5 118 Pr+5 119 } 



IV. RESULTS AND DISCUSSION 

To observe the physical situation of the problem 
under study, the velocity field, temperature field, 
concentration field, skin-friction, rate of heat transfer and 
rate of mass transfer are discussed by assigning numerical 
values to the parameters encountered into the corresponding 
equations. To be realistic, the values of Schmidt number 
(Sc) are chosen for hydrogen (Sc=0.22), helium (5c=0.30), 
water-vapor (5c=0.60), ammonia (5c=0.78), carbondioxeid 
(5c=0.96) at 25°C and one atmosphere pressure. The values 
of Prandtl number (Pr) are chosen for air (Pr=0.71), 
ammonia (Pr=0.90), carbondioxied (Pr=2.2), water 
(Pr=7.0). Grashof number for heat transfer is chosen to be 
Gr=10.0, 15.0, -10.0, -15.0 and modified Grashof number 
for mass transfer Gm= 15.0, 20.0, -15.0, -20.0. The values 
Gr>0 and Gm>0 correspond to cooling of the plate while the 
value Gr<0 and Gr<0 correspond to heating of the plate. 

The values of magnetic parameter (M=0.0, 0.5, 1.5, 
2.0), suction parameter (f w =3.0, 5.0, 7.0, 9.0), Soret number 
(5*0=0.0, 2.0, 4.0, 8.0) and heat source parameter (5=0.0, 4.0, 
8.0, 12.0) are chosen arbitrarily. 

The velocity profiles for different values of the 
above parameters are presented from Fig.-l to Fig.-ll. All 
the velocity profiles are given against rj. 

In Fig.-l the velocity distributions are shown for 
different values of magnetic parameter M in case of cooling 
of the plate (Gr>0). In this figure we observe that the 
velocity decreases with the increases of magnetic parameter. 

The reverse effect is observed in the Fig. -2 in case 
of heating of the plate (Gr<0). The velocity distributions 
are given for different values of heat source parameter S in 
case of cooling of the plate in the Fig. -3. From this figure 
we further observe that velocity decreases with the increases 
of heat source parameter. The opposite phenomenon is 
observed in case of heating of the plate in Fig. -4. In Fig.-5 
the velocity distributions are depicted for different values of 
suction parameter f w in case of cooling of the plate. Here we 
see that velocity decreases with the increases of suction 
parameter. The reverse effect is seen in case of heating of 
the plate in Fig. -6. 

The Fig. -7 represents the velocity distributions for 
different values of Schmidt number Sc in case of cooling of 
the plate. By this figure it is noticed that the velocity for 
ammonia (5c=0.78) is less than that of hydrogen (5c=0.22) 
decreases with the increases of Schmidt number. The 
reverse effect is obtained in case of heating of the plate in 
Fig. -8. The Fig. -9 exhibits the velocity distributions for 
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different values of Soret number S 0 in case of cooling of the 
plate. Through this figure we conclude that the velocity 
increases increasing values of Soret number. The opposite 
phenomenon is seen in case of heating of the plate from 
Fig.- 10. 

In Fig.- 11 the velocity distributions are presented 
for different values of M, S, Sc, Pr and f w in both case of 
cooling and heating of the plate where S 0 =0.0 and U 0 =1.0. 

It is noted that for externally cooled plate (a) an 
increase in M decreases the velocity field is observed by 
curves (/) and (ii); (b) an increase in f w decreases the 
velocity field observed by curves (i) and (Hi); (c) an increase 
in S decreases the velocity field observed by curves (Hi) and 
(iv); (d) an increase in Sc decreases the velocity field 
observed by curves (i) and (v); (e) an increase in Pr 
decreases the velocity field is seen in the curves (/) and (vi); 
(/) all these effects are observed in reverse order for 
externally heated plate. All the velocity profiles attain their 
peak near the surface of the plate and then decrease slowly 
along v-axis. i.e. far away from the plate. It is clear that all 
the profiles for velocity have the maximum values at rj=l. 

The temperature profiles are displayed from Fig.- 
12 to Fig.- 14 for different values of Pr, S and/ w against rj. 

We see in the Fig.- 12 the temperature distributions 
for different values of Prandtl number Pr. In this figure we 
conclude that the temperature decreases with the increases 
of Prandtl number. This is due to fact that there would be a 
decrease of thermal boundary layer thickness with the 
increase of prandtl number. In Fig.- 13 displays the 
temperature distributions for different values of heat source 
parameter S. In this figure it is clear that the temperature 
field increases with the increases of heat source parameter. 

Fig.- 14 is the temperature distributions for different 
values of Pr, S, f w for Gr=10.0, Gm=15.0, M=0.0, S 0 =0.0 
and Sc=0.22. In this figure it is obtained that an increases in 
Pr and S decreases the temperature field but an increases in 
of/ w increases the temperature field. 

The profiles for concentration are presented from 
Fig.- 15 to Fig.-17 for different values of So, Sc and/ w . 

The Fig.- 15 demonstrates the concentration 
distributions for different values of Soret number. It is 
noticed that the concentration profile increases with the 
increases of Soret number. Fig.-16 gives the concentration 
distributions for different values of suction parameter f w . We 
are also noticed that the effects of increasing values of f w 
decrease the concentration profile in the flow field. In Fig.- 
17 the concentration distributions are shown for different 
values of Sc, S 0 mdf w for Gr=10.0, Gm=15.0, M=0.0, S=2.0 
and Pr=0.71. In this figure it is seen that the concentration 
field increases with the increases of Schmidt number Sc. 

The concentration field also increases with the 
increases of Soret number S 0 . The concentration field 
decreases with the increases of suction parameter f m All the 
concentration profiles reach to its maximum values near the 
surface of the vertical plate i.e. rj=l and then decrease 
slowly far away the plate i. e. as //— k» 

The numerical values of skin-friction (x) at the 
plate due to variation in Grashof number (Gr), modified 
Grashof number (Gm), heat source parameter (S), Soret 
number (So), magnetic parameter (M), Schmit number (Sc), 
suction parameter (f w ), and Prandtl number (Pr) for 
externally cooled plate is given in Table- 1. It is observed 
that both the presence of S and So in the fluid flow decrease 



the skin-friction in comparison to their absence. The 
increase of M, Sc and f w decreases the skin-friction while an 
increase in Pr, Gr and Gm increase the skin-friction in the 
absent of S and S 0 . 

Table-2 represents the skin-friction for heating of 
the plate and in this table it is clear that all the reverses 
phenomena of the Table- 1 are happened. 

The numerical values of the rate of heat transfer in 
terms of Nusselt number (Nu) due to variation in Pandit 
number (Pr), heat source parameter (S) and suction 
parameter (f w ) is presented in the Table-3. It is clear that in 
the absence of heat source parameter the increase in Pr 
decreases the rate of heat transfer while an increase in f w 
increases the rate of heat transfer. In the presence of heat 
source parameter the increasing values of f w also increases 
the rate of heat transfer while an increase in Pr decreases the 
rate of heat transfer. It is also clear that the rate of heat 
transfer decreases in presence of S than the absent of S. 

Table-4 depicts the numerical values of the rate of 
mass transfer in terms of Sherwood number (Sh) due to 
variation in Schmit number (Sc), heat source parameter (S), 
Soret number (S 0 ), suction parameter (f w ) and Pandit number 
(Pr). An increase in Sc or Pr increases the rate of mass 
transfer in the presence of S and S 0 while an increase in f w 
decreases the rate of mass transfer. 



V. FIGURES AND TABLES 
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Fig.-l. Velocity profiles when Gr= 10.0, Gm= 15.0, 
5=0.0, S 0 =0.0, Sc=0.22, Pr=0.71,/ W =5.0 and U 0 =l.O 
against tj. 
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Fig. -2. Velocity profiles when Gr= -10.0, Gm= -15.0, Fig. -6. Velocity profiles when Gr= -10.0, 

5=0.0, S 0 =0.0, S*c=0.22, Pr=0.71,/ W =5.0 and U o =1.0 M=0.0, S 0 =0.0, S=0.0, 5*c=0.22, Pr=0.1l 

against rj against r\. 

-3=0.0 



Gm= -15.0, 
and U o=l.O 




12 3 4 5 6 

Fig. -3. Velocity profiles when Gr=10.0, Gm= 15.0, 
M=0.0, 5*0=0.0, S*c=0.22, Pr=0.71,/ W =5.0 and U 0 =l.O 
against rj. 

> 5=0.0 




Fig. -4. Velocity profiles when Gr=-10.0, Gm=-15.0, 
M=0.0, 5*0=0.0,^=5.0, 5*c=0.22, Pr=0.71 and U 0 =\.0 
against rj. 

.... ■ fc=3.0 




2 | (S S 10 

Fig. -5. Velocity profiles when Gr= 10.0, Gm= 15.0, 
M=0.0, S 0 =0.0, 5*=0.0, 5*c=0.22, Pr=0.71 and U 0 =l.O 
against rj. 





Fig. -7. Velocity profiles when Gr= 10.0, Gm= 15.0, 
M=0.0, S 0 =0.0, 5*=0.0, f w =3.0, Pr=0.71 and U 0 =l.O 
against r\. 



1 








i ■ ■ ■ ^0.22 




w-e-*.Si;=0.3Q 








5c=0.96 


Velocity profiles when Gr= 
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Fig.-8. 

M=0.0, 5*o=0.0, 5*=0.0,/ w =3.0, Pr=0.71 and C/^l.O 
against r\ 




2 4 0 

Fig. -9. Velocity profiles when Gr= 10.0, Gm= 15.0, 
M=0.0, 5*=0.0, 5*c =0.22, f w =5.0, Pr=0.71 and U 0 =l.O 
against rj. 
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Fig. -10. Velocity profiles when Gr= -10.0, Gm= -15.0, 
M=0.0, 5=0.0, / w =5.0 , Sc =0.22 



Ciwliii' &f flit plate (0=10.0. &n=ls.tt) M & Sc Pr £ 




Bering of tiepkte (&= -11.1 &i=-15.0) 



Fig.-l 1. Velocity profiles for different values of M, 5, 5c, 
Pr and f w when 5 0 =0.0 and t/ 0 =1.0 in both heating and 
cooling plate against rj. 




Fig. -12. Temperature profiles when Gr=10.0, Gm=15.0, 
M=0.5, 5=1.0, S 0 =4.0, 5c=0.22 and f w =3.0 against rj. 
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Fig.- 13. Temperature profiles when Gr= 10.0, Gm= 15.0, 
M=0.0, 5 0 =1.0, 5c=0.22, Pr= 0.71 and/ w =3.0 against rj. 

4 




Fig.- 14. Temperature profiles when Gr= 10.0, Gm= 15.0, 
M=0.0, 5 0 =0.0 and 5c=0.22 against rj. 




Fig. -15. Concentration profiles when Gr= 10.0, Gm= 15.0, 
5=0.0, M=0.5, 5c=0.22, Pr=0.71 and / w =5.0 against rj 




Fig. -17. Concentration profiles for different values of Sc, 
S 0 and/ w when Pr=0.71, M=0.0 and 5 =2.0 against rj. 



Table- 1: Numerical values of Skin-Friction (r) due to 
cooling of the plate. 



s. 

No 


Gr 


Gm 


S 


So 


M 


Sc 


f» 


Pr 


T 


1. 


10.0 


15.0 


0.0 


0.0 


0.5 


0.22 


5 


0.71 


3.9187 


2. 


15.0 


15.0 


0.0 


0.0 


0.5 


0.22 


5 


0.71 


4.3051 


3. 


10.0 


20.0 


0.0 


0.0 


0.5 


0.22 


5 


0.71 


4.6694 


4. 


10.0 


15.0 


4.0 


0.0 


0.5 


0.22 


5 


0.71 


3.7917 


5. 


10.0 


15.0 


0.0 


2.0 


0.5 


0.22 


5 


0.71 


4.6313 


6. 


10.0 


15.0 


0.0 


0.0 


2.0 


0.22 


5 


0.71 


3.2845 


7. 


10.0 


15.0 


0.0 


0.0 


0.5 


0.60 


5 


0.71 


2.5551 


8. 


10.0 


15.0 


0.0 


0.0 


0.5 


0.22 


7 


0.71 


2.5571 


9. 


10.0 


15.0 


0.0 


0.0 


0.5 


0.22 


5 


0.90 


4.1025 
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Table-2: Numerical values of Skin-Friction (r) due to 








cooling 


of the plate 






s. 


Gr 


Gm 


c 
.3 


So 


M 


be 


fw 


Pr t 


No. 


















1. 


-10 


-15 


0.0 


0.0 


0.5 


0.22 


5 


0.71 -2.131 


2. 


-15 


-15 


0.0 


0.0 


0.5 


0.22 


5 


0.71 -2.517 


3. 


-10 


-20 


0.0 


0.0 


0.5 


0.22 


5 


0.71 -2.881 


4. 


-10 


-15 


4.0 


0.0 


0.5 


0.22 


5 


0.71 -2.004 


5. 


-10 


-15 


0.0 


2.0 


0.5 


0.22 


5 


0.71 -2.843 


6. 


-10 


-15 


0.0 


0.0 


2.0 


0.22 


5 


0.71 -1.809 


7. 


-10 


-15 


0.0 


0.0 


0.5 


0.60 


5 


0.71 -0.767 


8. 


-10 


-15 


0.0 


0.0 


0.5 


0.22 


1 


0.71 -0.662 


9. 


-10 


-15 


0.0 


0.0 


0.5 


0.22 


5 


0.90 -2.314 


rable-3: Numerical values of the Rate of Heat Transfer (Nu 


S. No. 




Pr 




s 




fw 




Nu 


1. 




0.71 




0.0 




3 




-1.3753 


2. 




0.71 




4.0 




3 




-1.4886 


3. 




7.0 




0.0 




3 




-7.2539 


4. 




7.0 




4.0 




3 




-7.6531 


5. 




0.71 




0.0 




5 




-0.8068 


6. 




0.71 




4.0 




5 




-0.9232 



Table-4: Numerical values of the Rate of Mass Transfer (Sh) 



s. 

No. 


Pr 


S 


So 




Sc 


Sh 


1. 


0.71 


0.0 


0.0 


3 


0.22 


-0.5214 


2. 


0.71 


4.0 


2.0 


3 


0.22 


7.4912 


3. 


0.71 


0.0 


2.0 


3 


0.22 


6.5525 


4. 


0.71 


4.0 


2.0 


3 


0.60 


111.935 


5. 


0.71 


4.0 


2.0 


5 


0.22 


1.0846 


6. 


7.0 


4.0 


2.0 


3 


0.22 


662.888 


7. 


0.71 


8.0 


4.0 


3 


0.22 


25.3419 



VI. CONCLUSION 

In the present research work, combined heat and 
mass transfer effects on MHD free convection flow past a 
flat plate is presented. The results are given graphically to 
illustrate the variation of velocity, temperature and 
concentration with different parameters. Also the skin- 
friction, Nusselt number and Sherwood number are 
presented with numerical values for various parameters. In 
this study, the following conclusions are set out 

(1) In case of cooling of the plate (Gr>0) the velocity 
decreases with an increase of magnetic parameter, heat 
source parameter, suction parameter, Schmidt number 
and Prandtl number. On the other hand, it increases 
with an increase in Soret number. 

(2) In case of heating of the plate (Gr<0) the velocity 
increases with an increase in magnetic parameter, heat 
source parameter, suction parameter, Schmidt number 
and Prandtl number. On the other hand, it decreases 
with an increase in Soret number. 

(3) The temperature increase with increase of heat source 
parameter and suction parameter. And for increase of 
Prandtl number it is vice versa. 

(4) The concentration increases with an increase of Soret 
number and Schmidt number. Whereas it decreases with 
an increase of suction parameter. 

(5) In case of cooling of the plate (Gr>0) the increase of 
magnetic parameter, Schmidt number and suction 



parameter decreases the skin-friction. While an increase 
of Grashof number, Modified Grashof number and 
Prandtl number increases the skin-friction. In case of 
heating of the plate (Gr<0) all the effects are in reverse 
order. 

(6) The rate of heat transfer or Nusselt number (Nu) 
increase with the increase of suction parameter while an 
increase in Prandtl number it decreases. It also 
decreases in presence of S in comparison to absence of 
S. 

An increase in Schmidt number or Prandtl number 
increases the rate of mass transfer (Sherwood number) in the 
presence of heat source parameter and Soret number while 
an increase in suction parameter f w decreases the rate of 
mass transfer. We also see that both the increase in source 
parameter and Soret number increases the rate of mass 
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